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Abstract 

Both the eigenf unctions and the local density of states (DOS) near the pentagonal 
defects on the fullerenes was calculated analytically as a numerically. The results 
shows that the low-energy DOS has a cusp which drops to zero at the Fermi energy 
for any number of pentagons at the tip except three. For three pentagons, the nonzero 
DOS across the Fermi level is formed. Graphite is an example of a layered material 
that can be bent to form fullerenes which promise important applications in electronic 
nanodevices. The spheroidal geometry of a slightly elliptically deformed sphere was 
used as a possible approach to fullerenes. We assumed that for a small deformation the 
eccentricity of the spheroid is much more smaller then one. We are interested in the 
elliptically deformed fullerenes C70 as well as in Ceo and its spherical generalizations like 
big C240 and C540 molecules. The low-lying electronic levels are described by the Dirac 
equation in (2-1-1) dimensions. We show how a small deformation of spherical geometry 
evokes a shift of the electronic spectra compared to the sphere and both the electronic 
spectrum of spherical and the shift of spheroidal fullerenes were derived. In the next 
study the expanded field-theory model was proposed to study the electronic states near 
the Fermi energy in spheroidal fullerenes. The low energy electronic wave functions 
obey a two-dimensional Dirac equation on a spheroid with two kinds of gauge fluxes 
taken into account. The first one is so-called K spin flux which describes the exchange 
of two different Dirac spinors in the presence of a conical singularity. The second 
flux (included in a form of the Dirac monopole field) is a variant of the effective field 
approximation for elastic flow due to twelve disclination defects through the surface of 
a spheroid. We consider the case of a slightly elliptically deformed sphere which allows 
us to apply the perturbation scheme. We shown exactly how a small deformation of 
spherical fullerenes provokes an appearance of fine structure in the electronic energy 
spectrum as compared to the spherical case. In particular, two quasi-zero modes in 
addition to the true zero mode are predicted to emerge in spheroidal fullerenes. The 
effect of a weak uniform magnetic field on the electronic structure of slightly deformed 
fuUerene molecules was also studied. It was shown how the existing due to spheroidal 



2 



deformation fine structure of the electronic energy spectrum splitted in the presence of 
the magnetic field. We shown that the fine structure of the electronic energy spectrum 
is very sensitive to the orientation of the magnetic field. We found that the magnetic 
field pointed in the x direction does not influence the first electronic level whereas it 
causes a splitting of the second energy level. Exact analytical solutions for zero-energy 
modes are found. The HOMO (highest occupied molecular orbital)-LUMO (lowest 
unoccupied molecular orbital) gap was calculated within the continuum field-theory 
model of fuUerenes. 



Chapter 1 
Introduction 



The purpose of this chapter is to present some examples of using topology and geome- 
try in a study of a new interesting class of carbon materials as fuUerenes. FuUerenes are 
cage-like molecules of carbon atoms. The discovery of these molecules has attracted 
considerable attention of both experimentalists and theorists due to their unique phys- 
ical properties. An further interest to carbon nanoparticles originates from the fact 
that the geometry is accompanied by topological defects. The topologically nontriv- 
ial objects play the important role in various physically interesting systems. The 't 
Hooft-Polyakov monopole in the non-Abelian Higgs model, instantons in quantum 
chromodynamics, solitons in the Skyrme model, Nielsen- Olesen magnetic vortices in 
the Abelian Higgs model(see, e.g., [1]) are examples of such topologically nontrivial 
objects. Notice that the topological origin the elastic flux results in appearance of the 
disclination-induced Aharanov-Bohm like phase (see [2j). In condensed matter physics 
for instance, vortices in liquids and liquid crystals, solitons in low-dimensional systems 
are objects with nontrivial topology. 

The high flexibility of carbon allows producing variously shaped carbon nanostruc- 
tures: fuUerenes, nanotubes, nanohorns, cones, toroids, graphitic onions, etc. Some of 
them as e.g.nanotubes are composition of tubes and fuUerenes (as in Fig. II. II) .Historically, 
the fuUerenes Cgo (nicknamed also as Buckminsterfullerene or 'bucky ball') were first 
discovered in 1985 [3]. They are tiny molecular cages of carbon having 60 atoms and 
making up the mathematical shape called a truncated icosahedron (12 pentagons and 20 
hexagons). The amount of Ceo had been produced in small quantities until 1990 when 
an efficient method of production was achieved Since then, there were produced 
variously shaped fullerene molecules. These molecules are formed when 12 pentagons 
are introduced in the hexagonal lattice of graphene. Soon after the fuUerenes, the 
carbon nanotubes of different diameters and helicity [5] were produced. We can also 
image the construction of nanotube with fuUerenes at the end of nanotubes 

(see Fig. ll.2p . The mechanical, magnetic, and especially electronic properties of carbon 
nanotubes are found to be very specific (see, e.g., [B]). 

The fuUerenes are composed of hexagonal and pentagonal carbon rings. However, 
the structures having heptagonal rings are also possible. A junction which connect 
carbon nanotubes with different diameters through a region sandwiched by pentagon- 
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C7S nanotube 



Figure 1.1: The symmetric structure of fullerene Ceo, elliptically deformed molecule 
C78 and big deformed structure as carbon nanotube. 




Figure 1.2: The fullerene Cgo and spheroidal ones as a cap at the end of some nanotubes. 

heptagon pair has been observed in the transmission electron microscope [7J. A pair 
of five- and seven-membered rings have to be imposed to connect two different types 
of carbon nanotubes |8j. Theoretically the closed fullerenes and nanotubes exhibiting 
high topologies (from genus 5 to genus 21) were suggested in [9]. This follows from 
the known Euler's theorem that relates the number of vertices, edges and faces of an 
object. For the hexagonal carbon lattice it can be written in the form [0] 

...2^4 + ^5 - - 2^8... = - x)n^ = X = 12(1 - g), (1.1) 

where rix is the number of polygons having x sides, x is the Euler characteristic which 
is a geometrical invariant related to the topology of the structure, and g is the genus 
or a number of handles of an arrangement. According to (11.11) there is no contribu- 
tion to the Gaussian curvature for x = 6. This means that two-dimensional carbon 
lattice consisting only of hexagons is fiat. By their nature, pentagons (as well as other 
polygons with x 7^ 6) in a graphite sheet are topological defects (see Fig. 11.31) . In par- 
ticular, fivefold coordinated particles are orientational 60° disclination defects in the 



Figure 1.3: The pentagonal and heptagonal defects in the hexagonal network. 



sixfold coordinated triangular lattice. Description of the electronic structure requires 
formulating a theoretical model describing electrons on arbitrary curved surfaces with 
disclinations taken into account. The most important fact found in fiUl [TTl [T2] is that 
the electronic spectrum of a single graphite plane linearized around the corners of the 
hexagonal Brillouin zone coincides with that of the Dirac equation in (2+1) dimensions. 
This finding stimulated a formulation of some field-theory models for Dirac fermions on 
hexatic surfaces to describe electronic structure of variously shaped carbon materials: 
fuUerenes fl3[ [H] and nanotubes [15] . The field-theory models for Dirac fermions on a 
plane and on a sphere [16] were invoked to describe variously shaped carbon materials. 
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Chapter 2 

Electronic properties 



Carbon nanoparticles and their electronic properties have been a subject of an exten- 
sive study. Electronic states in nanotubes, fuUerenes, nanocones, nanohorns as well as 
in other carbon configurations are the subject of an increasing number of experimental 
and theoretical studies. It has been predicted and later observed in experiment that 
bending or stretching a nanotube change its band structure changing therefore the elec- 
trical properties: stretched nanotubes become either more or less conductive. Carbon 
nanotubes can be either a metal or semiconductor, depending on their diameters and 
helical arrangement [U] . This finding could allow to build nanotube-based transducers 
sensitive to tiny forces. 

The electronic properties depends on the topological defects. The peculiar elec- 
tronic states due to topological defects have been observed in different kinds of carbon 
nanoparticles by scanning tunneling microscopy (STM). For example, STM images 
with five-fold symmetry (due to pentagons in the hexagonal graphitic network) have 
been obtained in the Ceo fullerene molecule [H]. The pecuhar electronic properties 
at the ends of carbon nanotubes (which include several pentagons) have been probed 
experimentally in [121 ED]- Recently, the electronic structure of a single disclination 
has been revealed on an atomic scale by STM |21j . 

The problem of peculiar electronic states near the pentagons in curved graphite 
nanoparticles was the subject of intensive theoretical studies in fullerenes [HI [T3] . 
nanotubes [I5], nanohorns [22], and cones [23]. In particular, an analysis within the 
effective-mass theory has shown that a specific y/S x -\/3 superstructure induced by 
pentagon defects can appear in nanocones [21|. This prediction has been experimentally 
verified in [21J. A recent study [25j within both tight-binding and ab initio calculations 
shows a presence of sharp resonant states in the region close to the Fermi energy. The 
localized cap states in nanotubes have been recently studied in [26] . 

It is interesting to note that the problem of specific electronic states at the Fermi 
level due to disclinations is similar to that of the fermion zero modes for planar sys- 
tems in a magnetic field. Generally, zero modes for fermions in topologically nontrivial 
manifolds have been of current interest both in the field theory and condensed matter 
physics. As was revealed, they play a major role in getting some insight into un- 
derstanding anomalies [27] and charge fractionalization that results in unconventional 
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charge-spin relations (e.g. the paramagnetism of charged fermions) [2Sj with some 
important imphcations for physics of superfluid hehum (see, e.g., review [29]). 3-D 
space-time Dirac equation for massless fermions in the presence of the magnetic field 
was found to yield — 1 zero modes in the N- vortex background field [30]. As it was 
shown in [131, the problem of the local electronic structure of fuUerene is closely related 
to Jackiw's analysis [SO]- An importance of the fermion zero modes was also discussed 
in the context of the high-temperature chiral superconductors [21] ■ Among different 
theoretical approaches to describe graphene compositions, the continuum model play a 
special role. It allow us to describe some integral features of similar carbon structures. 
It should be noted that the formulation of any continuum model begins from descrip- 
tion of a graphite sheet (graphene). We briefly report now the effective mass theory for 
graphene [TH |T2]. A single layer of graphite forms a two-dimensional material. Two 
Bloch functions constructed from atomic orbitals for two inequivalent carbon atoms, 
A and B, in the unit cell provide the basis function for graphene. The Bloch basis 
functions $j [j = A, B) are constructed from atomic 2pz orbitals 

$A(^,r) ^^''^"/(^"'-^n), (2.1) 



^^(k^f) =J2^ e"^ ''^fif-d-r-), (2.2) 



d is a coordinate of B atom in the unit cell, A^ is a number of unit cell, fn is a unit cell 
coordinate, / is 2pz orbital, k is the wave vector. Now we want to solve the Schrodinger 
equation 

H^ = (2.3) 

where $ = ca^a + cb^b and H is the Hamiltonian of the system. We get the secular 
equations 

HaaCa + HabCb = ECa, 

HbaCa + HbbCb = ECb, (2.4) 

here 

ifi, = 1$,), i,j = AB. (2.5) 

In a tight-binding approximation it have a form 



Haa — Hbb — ^2 



H^^ = t (e*^"i + e^^"^ + e*^"^) , (2.6) 

where orbital energy of the 2pz level and are the coordinates of the three 

nearest-neighbor B atoms relative to an A atom 

t={f{r)\H\f{f-d)). (2.7) 



Solving the secular equation 

det {H -E) = 0, 



(2.8) 
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the eigenvalues E{k) are obtained in the form 



E{k) = ± t V 1 + 4 COS ^ COS + 4 cos2 (2.9) 

We have 

E{K) = E{K') = £2p, (2.10) 

and the K and K points are the corners of the Brillouin zone. In k.p theory we 
approximate the wave function at wave vector k = K + khj 

$(fc, f) = CA{k)e'''-'<l>A{K, r) + CB{k)e'''-'<l>B{K, r). (2.11) 

Inserting $ into Schrodinger equation, keeping terms of order k, we get the secular 
equation {h = c = 1 units are used) 

1_ f PAA PAB \ ( cAk) \ ^ / CAiR) \ ^2.12) 

m \ Pba Pee J \ cb[k-) J ^ ' \ cb[k) J 

where E{k) = E{k) - E{K) and 

p^- = J ^*{k,r)p^j{k,r)dfi,j = A,B. (2.13) 

It can be shown from group theoretical argument (and can be verify directly within 
the one-orbital tight-binding model) that the secular equation can be write in the 
form [HI HD] 

where vp = -^^-|^, a is the lattice constant and 

E{k)=±vf\k\. (2.15) 

The upper half of the energy dispersion curves describes the 7r*-energy anti-bonding 
band, and the lower half is the vr energy bonding band. Since there are two vr electrons 
per unit cell, the lower tt band is fully occupied. The points K,K create the Fermi 
surface of the graphene. When an external gauge fields are imposed, the translation 
symmetry can be broken and eigenfunction can no longer be labeled by k,. This requires 
a generalization of the trial wave function 

$(f) = J CA{K)e'^-^<l>A{K,r)dK + J CB{k)e'^-^<l>B{k,f)dK. (2.16) 

Inserting this function into Schrodinger equation leads to an equation [ini [32] 

{-^v,a.iV-^W)-E)(lf^)=0, (2.17) 
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here a — {ax,ay) are conventional Pauli matrices, W is the U(l) external gauge field 
and ^A{r): ^B{r) are envelope functions, smoothly varying in the scale of the lattice 
constant a, multiplying graphene Bloch function 

$(f) = ^^(f)$^(i^, f) + ^B{r)^B{K, r), (2.18) 

^!i{f)^ j Ci{k)e'^-^dn i^A,B. (2.19) 

Eq.(2.17) is algebraically identical to a two-dimensional Dirac equation, where com- 
ponents of spinor represent localization of electron on the graphene sublattice A and 
sublatticc B, respectively. We can get the similar equation for envelope functions 
near K point. Formally Eq.(2.17) can be obtained from Eq.(2.14) by exchanging 
Kk -i{dk - iWk), ca{b) *A(s)(r) and E{k) E. 



Chapter 3 
Theory 



Now we will formulate a theoretical model describing electrons on arbitrary curved 
surfaces with disclinations taken into account to investigate the electronic structure 
of fuUerenes. As it was mention in the previous section the electronic spectrum of 
a single graphite plane linearized around the corners of the hexagonal Brillouin zone 
coincides with that of the Dirac equation in (2+1) dimensions. This finding stimu- 
lated formulation of some field-theory models for Dirac fermions on hexatic surfaces 
to describe electronic structure of variously shaped carbon materials: fuUerenes [T3], 
nanotubes [l5l[26], and cones [33l [3^. 

The effective-mass theory for a graphene sheets can be generalized to graphene 
surfaces containing five carbon atom rings and creating geometrical structures as a 
cones or fuUerenes. In the continuum description pentagonal ring are described as 
a localized fictitious gauge fluxes. In such structures there is no globally coherent 
distinction between K and K points. 

Following the ideas above and describe fermions in a curved background, we need a 
set of orthonormal frames {cq} which yield the same metric, g^^y, related to each other 
by the local 5*0(2) rotation, 

eo.^e'^ = kiep, G 50(2). 

It then follows that g^^ = e'^e^dap where is the zweibein, with the orthonormal frame 
indices being a,/? = {1,2}, and coordinate indices = {1,2}. As usual, to ensure 
that physical observables are independent of a particular choice of the zweinbein fields, 
a local so(2) valued gauge field must be introduced. The gauge field of the local 
Lorentz group is known as the spin connection. For the theory to be self-consistent, 
the zweibein fields must be chosen to be covariantly constant |35j : 

V,e: := d,e: - F^.e^ + (^^)^e^ = 0, 

which determines the spin connection coefficients explicitly 

{u,r^ = e^D^r (3.1) 

Finally, the Dirac equation on a surface S in the presence of the U{1) external gauge 
field Wfj^ is written as |36j 
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(V;. - la^ - tW^) V^(r) = E^{r). (3.2) 



where V^i = + il^ with 



^M = ^^"/[7.,7/3], (3.3) 
being the spin connection term in the spinor representation and 

r = -la. = - ( I ) (3.4) 

2n 

a^dx^ = —T2I. (3.5) 

The matrix T2 acting in i^'-spin space. The energy E is measured relative to the Fermi 
energy, and the spinor tp = (\E'^\E'^\E'^ where '^{r) are envelope functions. For 

our purpose, we need incorporating both a disclination field and a nontrivial back- 
ground geometry. A possible description of disclinations on arbitrary two-dimensional 
elastic surfaces is offered by the gauge approach [3^. In accordance with the basic 
assumption of this approach, disclinations can be incorporated in the elasticity theory 
Lagrangian by introducing a compensating U{1) gauge fields W^. It is important that 
the gauge model admits exact vortex-like solutions for wedge disclinations (2^ thus 
representing a disclination as a vortex of elastic medium. The physical meaning of 
the gauge field is that the elastic flux due to rotational defect (that is directly con- 
nected with the Frank vector (see the next section)) is completely determined by the 
circulation of the PV^ field around the disclination line. In the gauge theory context, 
the disclination field can be straightforwardly incorporated in (13. 2p by the standard 
substitution d^ = d^ — iW^. 

Within the linear approximation to gauge theory of disclinations (which amounts 
to the conventional elasticity theory with linear defects), the basic field equation that 
describes the U{1) gauge field in a curved background is given by 

where covariant derivative := d^ + T^ involves the Levi-Civita (torsion- free, metric 
compatible) connection 

with g^k being the metric tensor on a Riemannian surface S with local coordinates 
x^ = {x'^,x'^). For a single disclination on an arbitrary elastic surface, a singular 
solution to (13. 6p is found to be ^7] 

W' = -ue''^DxG{x,y), (3.8) 

where 

D^D^G{x\ x^) = 2tt5\x\ x^)/,/g, (3.9) 
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with efj_k = ^/g^fJ.k being the fully antisymmetric tensor on S, ei2 = —£21 = 1. It 
should be mentioned that eqs. (13 .6113. 91) self-consistently describe a defect located on 
an arbitrary surface [37] . 

The general analytical solution to (13.21) is known only for chosen geometries. One 
of them is the cone ^331 El] • For the sphere which are of interest here, there were used 
some approximations. In particular, asymptotic solutions at small r (which allow us 
to study electronic states near the disclination line) were considered in [23] ■ For this 
reason, also the numerical calculations were performed in |38]. The results of both 
analytical and numerical studies of spherical geometry will be presented in the next 
section. 
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Chapter 4 
Spherical molecules 



According to the Euler's theorem, the fullerene molecule consists of twelve disclinations. 
Generally, it is difficult to take into account properly all the disclinations. There are 
two ways to simplify the problem. The ffist one is considering a situation near a single 
defect. This approximation we briefly sketch in this chapter. The second possibility 
is replacing the fields of twelve disclinations by the fields of the magnetic monopole. 
This approximation will be described in the next section. 



4.1 The Model 

We start the section with the computation of disclination fields. It will be used as a 
gauge fields in Dirac equation on the sphere. To describe a sphere, we employ the polar 
projective coordinates = r, = y?; < r < oo, < ip < 2tt with R being the radius 
of the sphere. In these coordinates, the metric tensor becomes 

grr = AR^/{R^ + r^, g^^ = AR^r^iR^ + , g,.^ = g^r = 0, (4.1) 

so that 



^:=^detM=4R'r/{R' + r')'. 
Nonvanishing connection coefficients (13.71) take the form 

2r R2 — 1 — 



and the general representation for the zweibeins is found to be 

e\ = e\ = 2R^cosip/{R^ + r^), e\ = -e^^ = -2R^ sm^/{R^ + r^), 
which in view of (13. ip gives 

u;l' = u;^' = 0, u^f = -i.f = 2r^/{R^ + r^)=:2u. (4.2) 
From the solutions (13.81) and (13.91) one can easily found 

G' = logr; Wr = Q, W^ = u, r ^ 0. 
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It describes locally a topological vortex on the Euclidean plane, which confirms the 
observation that disclinations can be viewed as vortices in elastic media. The elastic 
flux is characterized by the Frank vector uj, \uj\ = 27rz/ with u being the Frank index 
and the elastic flow through a surface on the sphere is given by the circular integral 

— / Wdf = u. 
27r7 

There are no restrictions on the value of the winding number u apart from v > —1 
for topological reasons. This means that the elastic flux is 'classical' in its origin; i.e., 
there is no quantization (in contrast to the magnetic vortex). If we take into account 
the symmetry group of the underlying crystal lattice, the possible values of v become 
'quantized' in accordance with the group structure (e.g., u = 1/6, 1/3, 1/2, ... for the 
hexagonal lattice). We define v = n/6, where n is the number of pentagons at the 
apices. It is interesting to note that in some physically interesting applications vortices 
with the fractional winding number have already been considered. A detailed theory 
of magnetic vortices on a sphere has been presented in [39] . 

The Dirac matrices in 2D space, can be chosen as the Pauli matrices: 7^ = — o"^, 7^ = 

and (13. 3p reduces to 

= itoa^ (4.3) 

From the assumptions above the Dirac operator T> := i'^'^e^{W ^ + iW^) on the two- 
sphere becomes 



V = V^ = — 

2i?2 



e-'^(-dr + + - 



(4.4) 



The generator of the local Lorentz transformations G SO {2) takes the form 
—id^p, whereas the generator of the Dirac spinor transformations p(A) is 

S12 = ^[71,72] = ^(T^- 

For massless fermions serves as a conjugation matrix, and the energy eigenmodes 
are symmetric with respect to E = {cr^ipE = V'-s)- The total angular momentum of 
the 2D Dirac system is therefore given by 

Lz = -id^ + ^c^^ 

which commutes with the operator (14.41) . Accordingly, the eigenfunctions are classified 
with respect to the eigenvalues of Jz = 3 + 1/2, 3 = 0, ±1, ±2, and are to be taken 
in the form 

'''-(„m2:"..)- (^-^i 

From the (14. 4p is follows that the spin connection term can be taken into account 
by redefining the wave function as 

^ = + r2, (4.6) 



4.2. 
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which reduces eigenvalue problem (13.21) to 



drU — — ^-u = Ev, 

r 

-drv - ii^^in^i) = Eu, (4.7) 



where E = 2R^E/{R^ + r^). 



4.2 Extended Electron States 

We will consider an approximate solution to (14.71) . Because we are interested in elec- 
tronic states near the disclination line, we can restrict our consideration to the case of 
small r. In this solution to (14. 7p (with (14. 6 p taken into account) is found to be 

where 77 = ±(j — u), fj = ±(j — u + 1), and A is a normalization factor. Therefore, 
there are two independent solutions with ri{f]) > and 77(17) < 0- Notice that respective 
signs ± in (14. 8 p correspond to states with E > and E < 0. As already noted, 
serves as the conjugation matrix for massless fermions, and the energy eigenmodes are 
symmetric with respect to -E = 0. One can therefore consider either case, for instance, 
^ > 0. 

The important restrictions come from the normalization condition 

J{\u\^ + \v\'^)^dx^dx'^ = 1. (4.9) 

From (14. Sp . it follows that ~ E. On the other hand, the integrand in (14. 9 p must be 
nonsingular at small Er. This imposes a restriction on possible values of j. Namely, 
for ri,fi > one obtains j — u > —1/2, and for ri,fi < one has j — u < —1/2. As is 
seen, possible values of j do not overlap at any u. 

In the vicinity of a pentagon, the electron wave function reads 



(4.10) 



In particular, in the leading order, one obtains ^ ~ \fE-i ^ ~ ^i/3^-i/6^ g^j^^j v]/ ~ 
^1/6^,-1/3 fQj, ^ _ 1/6, 1/3, respectively. Because the local density of states diverges 
as r — > 0, it is more appropriate to consider the total density of states on a patch 
< r < 5 for small 5, rather than the local quantities. For this, the electron density 
should be integrated over a small disk |r| < 5 (recall that r,Lp are stereographically 
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Figure 4.1: The solutions v'{x),u'{x) for different n. 



projected coordinates on the sphere). The result is 



D{E,5) oc 



' iE6)5, z/ = 0; 

u = 1/6,5/6; 
{ESyl^d, z/=l/3,2/3; 
5, u= 1/2; 



(4.11) 



For the defect free case {v = 0) we obtain the well-known behavior of the total density 
of electronic states (DOS) in the 6 disk given by D{E,6) ~ ES"^ (in accordance with 
the previous analysis [ID]). For u = 1/6,1/3,2/3,5/6, the low-energy total DOS has 
a cusp which drops to zero at the Fermi energy. Most intriguing is the case where 
u = 1/2 and a region of a nonzero DOS across the Fermi level is formed. This implies 
local metallization of graphite in the presence of 180° disclination. In the fuUerene 
molecule, however, there are twelve 60° disclinations, and therefore, the case u = 1/6 
is actually realized. 



4.3 Numerical Results 

The numerical calculations for the different type of carbon nanostructures as well as 
for the case of sphere were presented in paper [38]. As a starting point, the analytical 
asymptotic solutions found in the previous section are considered. The initial value 
of the parameter r is defined as r = 10"'^. It is worth noting that the choice of the 
boundary conditions does not influence the behavior of the calculated wave functions 
and only the starting point depends on it. A dimensionless substitution x = Er is used. 
The normalized numerical solutions to (14. 7p are given in Fig. 14.11 The parameters are 
chosen to be = 0.01 and R= 1. Notice that here we present the solutions for dotted 
values f'(= v) and u'{= u). The local DOS is shown schematically in Figs. 14.2 1 14.31 for 
different n. The Fig. 14. 31 describes also the dependence of the local DOS on a position of 
the maximum value of integrand in (14. 9p (which actually characterizes the numerically 
calculated localization point of an electron). Here and below 5 = 0.1. Notice that 
in fact the choice of the value of S does not influence the characteristic behavior of 
LDOS. As is seen, the DOS has a cusp which drops to zero at the Fermi energy. The 
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Figure 4.2: Schematic densities of states near the Fermi energy in the case of sphere. 



case n — 3 becomes distinguished. Let us emphasize once more that in the fuUerene 
molecule there are twelve 60° disclinations, so that the case n = 1 is actually reahzed. 
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Figure 4.3: 3d schematic plotting of the DOS near the Fermi energy for n = 0,1,2 
(going from the front side to the back side). 



Chapter 5 

Spheroidal geometry approach to 
fuUerenes 

Geometry, topological defects and the peculiarity of graphene lattice have a pronounced 
effect on the electronic structure of fullerene molecules. The most extensively studied 
Cqq molecule is an example of a spherical fullerene nicknamed a "soccer ball" ^0\ . 
Moreover, the spherical fuUerenes as Ceo are stable towards fragmentation than the 
other bigger fuUerenes [H]. The family of icosahedral spherical fuUerenes is described 
by the formula 20(n^ + n/ + /^) with integer n and /. Other fuUerenes are either slightly 
(as C70) or remarkably deformed and their general nickname is a "rugby ball". The 
electronic structure of C70 cluster of D^h geometry has been studied in [12] by using 
the local-density approximation in the density-functional theory. It was clearly shown 
that the spheroidal geometry is of a decisive importance for observed peculiarities in 
electronic states of the Cm. What is important, the calculated energy levels were found 
to be in good agreement with photoemission experiments. Spheroidal fuUerenes can be 
considered as a initially flat hexagonal network which has warped into closed monosur- 
face by twelve disclinations. We are interested here in the slightly elliptically deformed 
big fuUerenes like C240 and C540 molecules, see Figures [^TT| 15.21 It should be noted 
that the formulation of any continuum model begins from the description of a graphite 
sheet (graphene) (see, e.g., [101 ESI EH HSJ [121 Hi]). The models which involve the influ- 
ence of both geometry and topological defects on the electronic structure by boundary 
conditions were developed in Refs. [ISl SS] • A different variant of the continuum model 
within the effective-mass description for fullerene near one pentagonal defect was sug- 
gested in previous section. There are attempts to describe electronic structure of the 
graphene compositions by the Dirac-Weyl equations on the curved surfaces where the 
structure of the graphite with pentagonal rings taken into account is imposed with the 
help of the so-called K spin contribution W7\ . In the case of spherical geometry, 
this is realized by introducing an effective field due to magnetic monopole placed at 
the center of a sphere [13]. In this section we use a similar approach. The difference 
is that we consider also the elastic contribution. Indeed, pentagonal rings being the 
disclination defects are the sources of additional strains in the hexagonal lattice. More- 
over, the disclinations are topological defects and the elastic flow due to a disclination 
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Figure 5.2: Picture of molecular model of giant spheroidal fullerene C540. 

is determined by the topological Frank index. For this reason, this contribution ex- 
ists even within the so-called " inextensional" limit which is rather commonly used in 
description of graphene compositions (see, e.g., [IE])- 

Recently, in the framework of a continuum approach the exact analytical solution for 
the low energy electronic states in icosahedral spherical fuUerenes has been found [19] . 
The case of elliptically deformed fullerenes was studied in [HO] where some numerical 
results were presented. In this section, we suggest a similar to [IH] model with the 
Dirac monopole instead of 't Hooft-Polyakov monopole for describing elastic fields. We 
consider a slightly elliptically deformed sphere with the eccentricity of the spheroid 
e <C 1. In this case, by analogy with [51] we use the spherical representation for the 
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eigenstates, with the shght asphericity considered as a perturbation. This allows us to 
find explicitly the low-lying electronic spectrum for spheroidal fuUerenes. 

The problem of Zeeman splitting and Landau quantization of electrons on a sphere 
was studied in Ref. [52] • To this end, the Schrodinger equation for a free electron on 
the surface of a sphere in a uniform magnetic field was formulated and solved. Our 
studies also cover slightly elliptically deformed fuUerenes infiuent by the weak uniform 
external magnetic field pointed in the different directions. 



5.1 The model 

Let us start with introducing spheroidal coordinates and writing down the Dirac oper- 
ator for free massless fermions on the Riemannian spheroid S"^. The Euler's theorem 
for graphene requires the presence of twelve pentagons to get the closed molecule. In a 
spirit of a continuum description we will extend the Dirac operator by introducing the 
Dirac monopole field with charge G inside the spheroid to simulate the elastic vortices 
due to twelve pentagonal defects. The K spin fiux which describes the exchange of two 
different Dirac spinors in the presence of a conical singularity will also be included in 
a form of t'Hooft-Polyakov monopole with charge g. 

The Dirac equation on a surface S in the presence of the abelian magnetic monopole 
field and the external magnetic field is written as [HI HTl HQ] 

^7"e/[V^ - iW^ - - tA^]^ = E^. (5.1) 

We impose the K spin fiux by the gauge field. In Cartesian coordinates 

'^gyz -2gxz 

T2, tty = —— r-Ta, = 0. (5.2) 



r(r^ — z'^) ' r(r2 — z'^ 

Disclination field we choose in the form 



-Gv Gx 
= Wy = -^—-y = 0, (5.3) 

r[r + z) r[r + z) 

in the region Rj^f (northern hemisphere) and 

= Wy = = 0, (5.4) 

r(r — z) r(r — z) 

in the region Rs (southern hemisphere). Here r = y/ x'^ + y'^ + z'^ . 



5.2 Electronic states near the Fermi energy in weak 
magnetic field pointed in the z direction 

Firstly we will study the structure of electron levels near the Fermi energy infiuent 
by the weak magnetic field pointed in the z direction (see Fig. 15.31) so that A = 
B {y, —X, 0) /2. The elliptically deformed sphere or a spheroid 



24 CHAPTER 5. ELECTRONIC STRUCTURE OF SPHEROIDAL FULLERENES 




Figure 5.3: The schematic picture of the spheroidal fuUerene in a weak uniform mag- 
netic field pointed in the z direction. 

^ + 4 + 4 = 1' (5-5) 
a'^ & 

may be parameterized by two spherical angles = (f), = 9 that are related to the 
Cartesian coordinates X, y ^ z cis follows 

x = a sin 6^ cos 0; ?/ = a sin 6' sin 0; z = c cos 9. (5.6) 

The nonzero components of the metric tensor for spheroid are 

gtf,^ = sin^ 9\ gee = cos^ 9 ^ (? sin^ 9^ (5.7) 

where a,c>0, 0<^^<7r, O<0< 27r. Accordingly, orthonormal frame on spheroid 
is ^ ^ 

ei = : — - ds] 62 = , = de , (5.8) 

a sm^ Va2 cos^ ^ + sin^ ^ ^ ' 

and dual frame reads 



= a sin 9 d(j); = y cos"^ 9 + sin^ 9 d9 . (5.9) 
A general representation for zweibeins is found to be 



e'^^ = asin9; e^g = 0; e^^ = 0; e = y cos^ 9 + sin^ 9 . (5.10) 

Notice that e'^^ is the inverse of e"^. The Riemannian connection with respect to the 
orthonormal frame is written as 



de' = -uj\ Ae^= ^^"'^ === d0 A e^ (5.11) 
V cos^ y + c"^ sm 9 

de^ = -uj\ A e' = 0. (5.12) 

Here A denotes the exterior product and d is the exterior derivative. From Eqs. 15.111 
and 15.121 we get the Riemannian connection in the form 

^^2 = -^^1 = / . on : 9 ■ 2n ' ^^2 = ^01 = 0- (5.13) 

V cos^ 9 + sm^ 9 
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We assume that the eccentricity of the spheroid e = y |1 — (c/a)^| is small enough. In 
this case, one can write down c = a + 5a where 5 (| 5 1) is a small dimensionless 
parameter characterizing the spheroidal deformation from the sphere. So, we can follow 
the perturbation scheme using 5 as the perturbation parameter. Within the framework 
of the perturbation scheme the spin connection coefficients are written as 

u;i2 = -cjji ^ cos 0(1-5 sin^^), (5.14) 

where the terms to first order in 5 are taken into account. The inclusion of the K-spin 
connection can be performed by considering two Dirac spinors as the two components 
of an SU{2) color doublet ip = {■ipi'ipi)'^ (see Ref. jH] for details). In this case, the 
interaction with the color magnetic fields in a form of nonabelian magnetic fiux ('t 
Hooft'-Polyakov monopole) is responsible for the exchange of the two Dirac spinors 
which in turn corresponds to the interchange of Fermi points. As was shown in Ref. jl6j . 
the Dirac equation for ip can be reduced to two decoupled equations for ijj^ and ipi, 
which include an Abelian monopole of opposite charge, g = ±3/2. 

The elastic fiow through a surface due to a disclination has a vortex-like struc- 
ture [13] and can be described by Abelian gauge field W^. Similarly to Ref. |14j, we 
replace the fields of twelve disclinations by the effective field of the magnetic monopole 
of charge G located at the center of the spheroid. The circulation of this field is deter- 
mined by the Frank index, which is the topological characteristic of a disclination. Try- 
ing to avoid an extension of the group we will consider the case of the Dirac monopole. 
The another possibility is to introduce the non-Abelian 't Hooft-Polyakov monopole 
(see Ref. [IH])- As is known, the vector potential VF^ around the Dirac monopole have 
singularities. To escape introducing singularities in the coordinate system let us divide 
the spheroid (similarly as it was done for the sphere in [5l]) into two regions, Rn and 
Rs, and define a vector potential in R]y and a vector potential in Rs- Notice 
that below includes both gauge fields and defined in Eq. 15. 1[ 

One has 

Rn : < < ^ + A; < < 27r, 
Rs : ^ - A < < tt; < < 27r, 

Rns ■ ^ - A < e < ^ + A; < (fx 2tt, (overlap) (5.15) 

where A is chosen from the interval 0<A<7r/2. In spheroidal coordinates (15.71) the 
only nonzero components of W^^ are found to be 

cos 6(1 + 6 sin^ 0) + G{1 - cos 6) - 6G sin^ 6 cos 6, (5.16) 

^g cos 6(1 + 6 sin^ 6) - G{1 + cos 6) - 6G sin^ 6 cos 6, (5.17) 

where the terms of an order of 5^ and higher are dropped. In the overlapping region 
the potentials and are connected by the gauge transformation 

= W^ + tSNsd^S^l (5.18) 
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where 

JNS — e 

is the phase factor. The wavefunctions in the overlap R^s are connected by 

i/^N = Sns^^s, (5.19) 

where ipN and ips are the spinors in R^ and Rs, respectively. Since Sns must be single- 
valued [55J, 2G takes integer values. Notice that the total flux due to Dirac monopole 
in Egs. 15.161 and 15.171 is equal to AirG. What is important, there is no contribution 
from the terms with 6. In our case, the total flux describes a sum of elastic fluxes due 
to twelve disclinations, so that the total flux (the modulus of the total Frank vector) 
is equal to 12 x tt/S = An. Therefore, one obtains G = 1. 

Let us consider the region Rjsf. In spheroidal coordinates, the only nonzero compo- 
nent of in region Rjy is found to be 

W^^g cos ^(1 + 6 sin^ 6) + ^(1 - cos 6) - 6G sin^ 6 cos 9, (5.20) 

and the external magnetic field pointed in the z-direction reads 

A^ = -^Ba^sin^e. (5.21) 

By using the substitution fl5.53p we obtain the Dirac equation for functions Uj and f j 
in the form 



where 



Pi = -— sin0(j - 2mcos^) -71— sin^0. (5.23) 
a 2 

A convenient way to study the eigenvalue problem is to square Eq. (15.221) . For this 
purpose, let us write the Dirac operator in Eq. fl5.22p as P = Vq^^Vx. One can easily 
obtain that 



-r, 1 cos6'^ 1 1 \ (j — m cos 6) ^ B , . „ 



m—j cos 9 B'^a^ sm^ 9 

+^3 2n + 1 • 5.24 

sm 9 4 

Notice that to first order in 6 the square of the Dirac operator is written as = 
("Dq + 6r), where T = {'Dq'Di + ViVq). In an explicit form 

a^T = 2j^ + jx ((J3 — 6m) + Am{m — (Ta)^^ + 2mcr3 + 3i?a^(l — x^)x{a'i/2 — m) 
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+ 2Bja\l - x^) + B^a^{l - x^)72, (5.25) 

where the appropriate substitution x = cos 9 is used. The equation 'D'^ip = E'^ip takes 
the form 

where A = aE, V{x) = j + {o'^ — 2m)x/2. Since we consider the case of a weak magnetic 
field the terms with and 6B in Eq. fl5.26p can be neglected. As a result, the energy 
spectrum for spheroidal fuUerenes is found to be 



(Ajj2 = + |^-| + 1/2)2 _ ^2 ^ ^^2 ^ ^ ^(2j2 + + A^) (5.27) 

j(m2 - 1/4) 



where 



p{p+l) 

and 

^ j^(6m^ - 1/2) ^ jW(4m^-3) ^ ^8) 

p(p + 1) p'^{p + ly 

^ (4m(m - l)F„(2a, 2/?) - X) T J„(2a, 2/j) + (4m(m + 1)F„(2^, 2u) + X) J„(2/i, 2u) 

T J„(2«,2/?) + J„(2/i,2i/) 

(5.29) 

p(p + 1) p^(p + 1)^ 

(n + l)(n + 2a + l)(n + 2/3 + l)(n + 2a + 2/? + l) 
= (2p + l)(p+l)2(2^ + 3) ^'-'^^ 

^n(n + 2a){n + 2/5)(n + 2a + 2(3) 



(2p- l)j92(2p+ 1) 
J„(2a, 2/5) = r(n + 2a + l)r(n + 2/5+1), 
Here p = n + P + a and 

T = I + (5.32) 



When S = 0, B = 0, one has the case of a sphere with magnetic monopole inside. 
In this case, the exact solution is known (see, e.g., [49] ) 



d,il-x')d^- 



2 (j - mxf - ja^x + 4 + cxsm 



l-x2 



<(x) ) ^^0" 4^ 1, <(x) 
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u 



jn 



with the energy spectrum 



a(l -a;)"(l + x)^p2a,2/3( 
a(l-a;)^(l + a;)'^p2/.,2^l 



x) 



{Xir = {n+\j\ + 1/2)'- 



(5.33) 



(5.34) 



Here 



a 



1 

2 
1 

2 



J — m 



,/3 



j — m + 



V 



j + m + 



j + m — 



(5.35) 



p2a,2/3^^^ and P^'^'^'^(x) are Jacobi polynomials of n-th order, and Cu and C^, are the 
normalization factors. These states are degenerate. For example, the degeneracy of 
the zero mode is equal to six. Therefore, we have to use the perturbation scheme for 
the degenerate energy levels [56] . 

Finally, in the linear in 6 approximation, the low energy electronic spectrum of 
spheroidal fullerenes takes the form 



Ejn 



SB, 



rpv I 7710, 



with 



(5.36) 



± 



E^ 



SB, 
jn 



Ba\j + Ajn) 



{2^ + n){2r] + n) 
6{2f + A' + A') 



2E^ 



(5.37) 



where ^ = (i^) and t] = (3 (a) for j > (j < 0), respectively. The parameters a, /3, 
/i, v are defined above. 

Here we came back to the energy variable E = A/a (in units of hVp/a where Vp 
is the Fermi velocity). Notice that the non-diagonal matrix element of perturbation 
{e^^\{il:jn\T\ij-jn)\e-^^) turns out to be zero. Therefore, the states with opposite j do 
not mix and j remains a good quantum number as would be expected. As is seen from 
Eq. (I5.36p . for 5 = both eigenstates and eingenvalues are the same as for a sphere 
(cf. Ref. [in])- At the same time, Eqs. fl5.36p and (15.371) show that the spheroidal 
deformation gives rise to an appearance of fine structure in the energy spectrum. The 
difference in energy between sublevels is found to be linear in 6 which resembles the 
Zeeman effect where the splitting energy is linear in magnetic field. In addition we 
found a further splitting of the states with opposite j. 

As an example. Table [SH] shows all contributions (in compliance with Eq. (I5.36p ) to 
the first and second energy levels for YO-C240 (YO means a structure given in [571158]). 
Schematically, the structure of the first energy level is shown in Fig. 15.41 As is clearly 
seen, the first (double degenerate) level becomes shifted due to spheroidal deformation 
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YO - C240 J ^°„(el^) El^{meV){B = {}) Eg- (me\/)(Ea' = 0.1) 



72 = 0, ?7i = 


1/2 


1 


1.094 


10.5 


37.5 






1 

- 1 


1 .uy^ 


1 n ^ 


-lu.o 


n = 1, m = 


1/2 


1 


1.89 


8.8 


24.8 






-1 


1.89 


8.8 


-7.2 


n = 0, m = 


1/2 


2 


1.89 


28.4 


60.4 






-2 


1.89 


28.4 


-3.6 


n = 0, m = 


-5/2 


3 


1.89 


3 


27 






-3 


1.89 


3 


-21 



Table 5.1: The structure of the first and higher energy levels for YO-C240 fuUerene in uniform 
magnetic field. The hopping integral and other parameters are taken to be t = 2.5 eV and Vp — 3tb/2'h, 
b = I.45A, R = 7.O3A, SD = O.I7A, S = 0.024. b is average bond length, R [R =^ a) is average radius, 
SD is standard deviation from a perfect sphere (see Refs. [STIEH]), so that 6 = SD/ R. 

and splitted due to nontrivial topology. The uniform magnetic field provides the well- 
known Zeeman splitting. The difference between topological and Zeeman splitting is 
clearly seen. In the second case, the splitted levels are shifted in opposite directions 
while for topological splitting the shift is always positive. The more delicate is the 
structure of the second level which is schematically presented in Fig. 15.51 In this case, 
the initial (for 5 = 0) degeneracy of E^^ is equal to six. The spheroidal deformation 
provokes an appearance of three shifted double degenerate levels (fine structure) which, 
in turn, are splitted due to the presence of topological defects. The magnetic field is 
responsible to Zeeman splitting. 




Figure 5.4: The schematic picture of the first positive electronic level E^-^ for spheroidal 
fuUerenes in a weak uniform magnetic field pointed in the z direction. 

We admit that the obtained values of the splitting energies can deviate from the 
estimations within some more precise microscopic lattice models and density-functional 
methods (see, e.g. [59]). Therefore, we focused mostly on the very existence of physically 
interesting effects. For this reason, the values in both Table and the schematic pictures 
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Figure 5.5: The schematic picture of the second positive electronic level Ej^^ of the 
spheroidal fuUerenes in a weak uniform magnetic field pointed in the z direction. 

are presented with accuracy at about one percent of E^^. 
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5.3 Electronic states near the Fermi energy in weak 
magnetic field pointed in the x direction 

In the second case we introduce that the uniform external magnetic field B is chosen 
to be pointed in the x direction (see 15.61) so that A = B {0, —z,y) /2. 




Figure 5.6: The schematic picture of the spheroidal fullerene in a weak uniform mag- 
netic field pointed in the x direction. 

The nonzero component of Wfj_ in region are the same as in previous case but 
the external magnetic field are changed 

= — -Bac sin 9 cos 9 cos (j), (5.38) 

Ag = —-BacsiiKf), (5.39) 

where c = a + 6a and only the terms to the first order in 6 are taken into account. 
Following the calculations above with neglecting the terms B^ and 6B we obtain the 
Dirac equation for functions Uj and Vj in the form 

^ 1 / o , cot .1 ^ 2 ■ i\ "^2 /. a ^ ID 2 



—i(Ji — (dg-\ \-i-Ba sin0) H j— mcos^H — Ba sin ^ cos ^ cos I + 51^1 x 

^ a 2 2 asm9 \ 2 J J 

X ( "'III )=e( ^'f2 ) . (5.40) 



where 



Vi = — — sin^ (j — 2m cos ^) — 71 — sin^ 9cos9cos(j). (5-41) 
a 2 



We find the square of nonperturbative part of Dirac operator above in the form 
f)2 _ 1 /^;)2 , cos^ 1 1 \ {j-mcos9f 



a 



2 I ^ sin^ ' 4 AsinH a'^sinH 



+a, '^ ^^.Tf +BV{9,<j>), (5.42) 
sm 9 
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where the terms with and 6B was neglecting and 

T^/zi ,x ■/ ,o sin</), ^ cos0(l + sin^ 6') ..^ 1 . , 

V{9,(f)) = -i{cos(f)d^ —) cot 6 -as . i{de + - cot 9) sm (f) 

Zi sill. \j Zi 

m cos^ Q cos (/) / . X 

(5.43) 



sin ^ 

We denote V^^ = {4'j^\V{9,(j))\ipj^^), where ipj^ are the eigenfunctions of Vq. In the 

case of first energy mode the terms V}{^_i, Viji , Vil\, Vli^i are zero and therefore the 
weak magnetic field pointed in the x direction do not influent the first energy level. 

1 /2 1 /2 

For the second energy mode only nondiagonal terms V_{ _2, 2 were find different 
from zero. The following numerical values of the nondiagonal terms VJ^, was found: 

Vli^_2 = -2.1, V^i^ = 2.1. We use the wave functions for calculations of nondiagonal 
terms in the case of second energy level in the form: 




^-2(^,0) = ^^\/ 77 . 9 ' (5-44) 




,V^(..,).£J,/|f 'Vi(-2;^+l) |. (5,45) 



The low energy electronic spectrum of spheroidal fullerenes in this case takes the 
form 

E^^ = El + E'''% (5.46) 



where 



|2 



E'^^ = ^ ^ ^ • (5-47) 

jn 

We change 1, 2 to —1, —2 in the expression above to get the another energy levels in 
the degenerate energy state which is considered. Where the Tu are the diagonal matrix 
element of perturbation of Eq. 15.251 with 5 = 0. Table 15.21 shows all contributions to 
the first and second energy levels for YO-C240 fullerenes influent by the horizontal weak 
external magnetic field. Figures 15. 7^ 15.81 present first and second positive energy levels 
of spheroidal fullerenes in weak magnetic field pointed in the x direction in comparison 
with Figures 15.41 15.51 As is seen, there is a marked difference between the behavior of 
the first and second energy levels in magnetic field. Indeed, in both cases the energy 
levels become shifted due to spheroidal deformation. However, the uniform magnetic 
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field does not influence the first energy level. The splitting takes place only for the 
second level. We can conclude that there is a possibility to change the structure of the 
electronic levels in spheroidal fuUerenes by altering the direction of magnetic field. It 
would be interesting to test this prediction in experiment. The values in both Table 
and the schematic pictures are presented the same way as in the previous case with 
accuracy at about one percent of Ej^. 



YO 


— C240 




j 


EUeV) 


El{meV){B = 0) 


El^'^{meV){Ba'^ = 0.1) 


n = 


0, m = 


1/2 


1 
-1 


1.094 
1.094 


10.5 
10.5 


10.5 
10.5 


n = 


0, m = 


-5/2 


3 
-3 


1.89 
1.89 


3 
3 


3 
3 


n = 
n = 


0, m = 
1 


1/2 


2 
1 


1.89 
1.89 


28.4 
8.8 


53/-16/ 


n = 
n = 


1, m = 



1/2 


-1 
-2 


1.89 
1.89 


8.8 
28.4 


53/-16/ 



Table 5.2: The structure of the first and second energy levels for YO-C240 fuUerene in uniform 
magnetic field. The hopping integral and other parameters are taken to be t = 2.5 eV and Vp = 3i5/2?i, 
h = I.45A, R = 7.O3A, SD = O.I7A, 5 = 0.024. b is average bond length, R {R = a) is average radius, 
SD is standard deviation from a perfect sphere (see Refs. [S7l[S5j), so that 5 — SD/R. 



Figure 5.7: The schematic picture of the first positive electronic level Ej^ for spheroidal 
fuUerenes in a weak uniform magnetic field pointed in the x direction. 
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Figure 5.8: The schematic picture of the second positive electronic level E^^ with a 



mixture of quantum numbers j = 1,2 and j = —1, —2 for spheroidal fuUerenes in a 
weak uniform magnetic field pointed in the x direction. 

5.4 Zero-energy mode 

Now we will analyze the special case of an electron state at the Fermi level (so-called 
zero-energy mode), in this case the right side of equation 15.11 is equal to zero. We will 
use the projection coordinates in the form 

2i?V 2R^r R^-r^ , , 

where R and c are the spheroidal axles. The Riemannian connection with respect to 
the orthonormal frame is written as 



where 



-"«="« = 7(F^=^r=5;5^ ^;.=^?, = o. (5,49) 

The only nonzero component of the gauge field in region R^ for spheroidal 
fullerenes reads 

W^ = G + m^, (5.50) 

Notice that the monopole field in Eq. (15.501) consists of two parts. The first one 
comes from the K-spin connection term and implies the charge g = ±3/2 while the 
second one is due to elastic flow through a surface. This contribution is topological in 
its origin and characterized by charge G. For total elastic flux from twelve pentagonal 
defects one has G = 1. The parameter m = g — G is introduced in Eq. (15.501) . The uni- 
form external magnetic field B for simplicity is chosen to be pointed in the z direction. 
In projection coordinates the only nonzero component of is written as 
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Let us study Eq. fl5.1l) for the electronic states at the Fermi energy {E = 0). The Dirac 
matrices can be chosen to be the Pauh matrices, 71 = —(72,72 = — ci. By using the 
substitution 



E 



Vj{r) 



.J 



0,±1,±2,... 



we obtain 



R' 



2r(r2 + i?2) 
r2-/?2 



) vir) = 0, 



-r^' -"^X- - -^^^^ - 2r(r2 + R^) ■ 



uir) 



0, 



(5.53) 



(5.54) 



(5.55) 



where K = ^ (i?2 _ 7-2^ _|_ 4c2j.2_ assume that the eccentricity of the spheroid is 
small enough. In this case, one can write down c = R + 6R and 5 (| 5 |^ 1) is a 
small dimensionless parameter characterizing the spheroidal deformation. Therefore, 
one can follow the perturbation scheme using S as the perturbation parameter. To the 
leading in 6 approximation Eqs. f l5.54p and (15.551) are written as 



drVj{r) 



{R^ + r 



2\2 



r[(i?2 + r2)2-4/?2r25] 
2BR^r 



[i — m{l + 5) 



+ 



r-2 + R2 
R^ 



(r2 

+ mo 



R 



2\3 



r2 + R'^Y 



drUj{r) 



[(i?2 + r2)2 - 4i?2r25] 2r(i?2 + r-2 
{R' + r^f 



)vj{r), 



(5.56) 



'r[(i?2 + r2)2-4i?2r25] 
2BR'^r 



.1 A^'-R'fi 



r 



r2 + i?2 

2 j^2 



[(i?2 + r2)2-4i?2r25] 2r(i?2 + r2) 
The exact solution to Eqs. fl5.56p and f l5.57p is found to be 



(5.57) 



Vj{x) 



(a; + i)(i/2-m) /x + l-25- 2^5(J 



Uj{x) = (x + 1) 



x+l-26 + 2J6{6-l) 



((x + 1)^ -45x)"', (5.58) 



{m+l/2)^i{i+m-l/2) 



a; + 1 - 25 - 2 J5(5 - 1 



x + l-26 + 2J6{6-l 



[x + 1)2 - 46x 



(5.59) 

where x = r^/R^, a = {j/2)^5/{5 - 1) + BR^/A^5{5 - 1). The function u^- can be 
normalized if the condition m— 1/2 < j < 1/2— m is fulfilled. In turn, the normalization 
condition for vj has the form — 1/2— m < j < m+1/2. Finally, there are five normalized 
solutions for uj with j = 0, ±1, ±2 when m = —5/2, and one solution for vj with j = 
when m = 1/2. In the limit 5 — >^ we arrive at the solution for spherical fuUerenes 

y .(r) = ^{-l/2-J-™)(/?2 ^ ^2yi/2+m)g Jf^^ 



Uj{r) 



r 



(-1/2+i+m) 



(i?2 + r 



2\(l/2-m) 



(5.60) 
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Since normalization conditions do not depend on the parameter 6 they are the same 
as for spherical case. Thus the perturbation does not change the number of zero 
modes. Notice that Eq. f l5.60p can be obtained from Eqs. (15.581) and fl5.59p by using 
lim [1 + (x/a)]" = e^,a — > oo. Because of the different behaviour of the zero-energy 
modes there is also possibility to recognize the zero-eigenvalue states in experiment. 



5.5 The SU{2) algebra 



The angular-momentum operators for Dirac operator on the sphere S"^ with charge G 
and total magnetic monopole m looks as follows 

L, = -^{^^T^G), (5.61) 

f • • ,o ■ , COS 61/^ . "m \ cos(f) 

L^ = ismcpde + zcos0^— - =F «G - i + cr^ (5.62) 

sm fc' V cos 6 / z sm 

f • ,n • • ,cos6'/ _ . m \ sin0 

Ly = -I cos (pde + 2sm0^— [dsTiC - i + cr, . ^ , (5.63) 

sm V cos d J Ismu 

where — (-I-) sign correspond to north (south) hemisphere. These operators satisfy the 
standard commutations relations of the SU{2) algebra: 

EijkLjLk = iLi. (5.64) 

The square of Dirac operator and L^ may be diagonalized simultaneously 

Vl{B = 0) = P + ^-m^. (5.65) 

We can see that the operator V{d, (p) in the case when the magnetic field is pointed in 
the z direction can be expressed in the form 

V{e,<j))=L,+ (^^-rnjz. (5.66) 

In the case when the magnetic field is pointed in the x direction we have 

l^(^,</.) = -4-(y-m)x, (5.67) 

where x,z are Cartesian coordinate 

X = sin 6^ cos 2; = cos6'. (5.68) 

The square of Dirac operator and operator may be diagonalized simultaneously and 
their eigenvalues are interrelated 

{E, n\L^\E, n) = l{l + l)=E^-^ = {n+ \j\){n + \j\ + 1), (5.69) 



5.5. 
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where 

l = n+\j\. (5.70) 

We can use instead of n the angular momentum I = n + \ We can introduce the 
eigensteite of T>o{B = 0) in the form 



2TTn ^ {i~3)\ 

isgn{E){l - a;)-i(j'-™+i/2)(i + a;)-|(j+™"V2) _ ^yi-m+m^^ ^ 3,)(/+m-i/2) j 

5.71) 



where 



n = 2Vr(/ - m + 1/2)2 + r(/ + m + 1/2)2, 



and r are the Gamma functions (see Ref. [60]). When we act by L_ and L+ on 
find that 



where 



L+^^;,, = ^(/ + j + l)(Z-j)^;,,+i, 



-e ^ 5e - [d^TT'G - i- 



5.72) 
^ we 

5.73) 
5.74) 

5.75) 

„ I „ , . 5.76) 

sin 6 \ ' cos 6 J ' 2 sin ^ ^ ^ 

Now we would like to transform our formulae to the cartesian coordinates. The corre- 
sponding transformation rules for spinors are [60] 

mc = V^^, (5.77) 

and cartesian realization of operator L is 

Lc = V^LV, (5.78) 

with l^-matrices given by 



sin^^ \ cos 9 J 2sin^ 



V = exp {^^0^ exp (^^^ 

For example we present the Cartesian realization of operator 

r9 \ ■ ■ Aa : ■ fa ■ m \ ax 

[Lx)c = ism.<pde + icos<p^— =f zG - i + — - 

smfc* V cos 6 J 2 

that may be expressed in the form (m = G = 0) 



Ox 



{Lx)c = ypz - zpy + —. 



5.79) 



5.80) 



5.81) 
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Chapter 6 
Conclusion 



Electronic properties of fullerenes molecules have been discussed from a theoretical 
point of view. The effective-mass approximation is particularly suitable for under- 
standing global and essential features. In this scheme, the motion of electrons in 
fullerenes is described by Dirac-Weyl's equation. The geometry and topology is found 
to influence the main physical characteristics of graphite nanoparticles, first of all, their 
electronic properties. The topological defects (disclinations) appear as generic defects 
in closed carbon structures. 

We have considered the electronic states of spheroidal fullerenes provided the spheroidal 
deformation from the sphere is small enough. In this case, the spherical representation 
is used for describing the eigenstates of the Dirac equation, with the slight asphericity 
considered as a perturbation. The using of the perturbation scheme allows us to find 
the exact analytical solution of the problem. In particular, the energy spectrum of 
spheroidal fullerenes is found to possess the fine structure in comparison with the case 
of the spherical fullerenes. We have shown that this structure is weakly pronounced 
and entirely dictated by the topological defects, that is it has a topological origin. We 
found three twofold degenerate modes near the Fermi level with one of them being the 
true zero mode therefore our finding confirms the results of [42j that there can exist 
only up to twofold degenerate states in the C70. 

Notice that for spherical fullerenes {6 = 0) our results agree with those found in [19]. 
It is interesting that the predictions of the continuum model for spherical fullerenes are 
in qualitative agreement with tight- binding calculations [611 [621 [631 [M] . particular, 
the energy gap between the highest-occupied and lowest-unoccupied energy levels be- 
comes more narrow as the size of fullerenes becomes larger. It is important to keep 
in mind, however, that the continuum model itself is correct for the description of the 
low-lying electronic states. In addition, the validity of the effective field approximation 
for the description of big fullerenes is also not clear yet. Actually, this approximation 
allows us to take into account the isotropic part of long-range defect fields. For bigger 
fullerenes, one has to consider the anisotropic part of the long-range fields, the influ- 
ence of the short-range fields due to single disclinations as well as the multiple-shell 
structure. 

We have focused also on the structure of low energy electronic states of spheroidal 
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fuUerenes in the weak uniform magnetic field. For the states at the Fermi level, we 
found an exact solution for the wave functions. It is shown that the external mag- 
netic field modifies the density of electronic states and does not change the number 
of zero modes. For non-zero energy modes, electronic states near the Fermi energy of 
spheroidal fullerene are found to be splitted in the presence of a weak uniform mag- 
netic field. The case of the x-directed magnetic field was considered and compared 
with the case of the z-th direction. The z axis is defined as the rotational axis of the 
spheroid with maximal symmetry. The most important finding is that the splitting of 
the electronic levels depends on the direction of the magnetic field. Our consideration 
was based on the using of the eigenf unctions of the Dirac operator on the spheroid, 
which are also the eigenfunction of L^. Let us discuss this important point in more 
detail. 

In the case of a sphere there is no preferable direction in the absence of the magnetic 
field. The magnetic field sets a vector, so that the z-axis can be oriented along the field. 
In this case, one has to use such eigenfunctions of the Dirac operator which are also 
the eigenfunctions of L^. For the x-directed magnetic field, the eigenfunction of both 
the Dirac operator and L^ must be used. Evidently, the same results will be obtained 
in both cases. 

The situation differs markedly for a spheroid. The spheroidal symmetry itself as- 
sumes the preferential direction which can be chosen as the z-axis. In other words, the 
external magnetic field does not define the preferable orientation. The symmetry is 
already broken and, as a result, the case of the magnetic field pointed in the x direction 
differs from the case of the ^-directed field. For instance, there are no eigenfunctions 
which would be simultaneously the eigenfunctions of both the Dirac operator on the 
spheroid and L^. For this reason, the structure of the electronic levels is found to 
crucially depend on the direction of the external magnetic field (Figures 5.3-5.8). 

It should be mentioned that the zero-energy states in our model corresponds to the 
HOMO (highest occupied molecular orbital) states in the calculations based on the 
local-density approximation in the density-functional theory (see, e.g. [l2j). In partic- 
ular, the HOMO-LUMO energy gap is found to be about 1.1 eV for YO-C240 fullerene 
within our model (here LUMO means the lowest unoccupied molecular orbital) [651 ES] . 

The very big fullerenes like Cggo and C1500 become more deformed, faceted and 
can no longer form a free-electron model like the electronic shell [HZ], which was the 
assumption for this model. For these structures the deviation from the sphericity 
is larger when the pentagon defects are localized at the opposite poles. In the case 
when the poles are far away from each other we obtain the structure of nanotubes, 
and for the exact description some new model related to that proposed here should 
be used. In our opinion, our predictions are quite general for the fullerene family 
and are of interest for experimental studies. Finally, we think that the spheroidal 
geometry approach especially our established generic models could also be related to 
other physical problems with slightly deformed spherical structures that are common 
in the nature e.g. in astrophysics for modeling pulsars or kvasars. 
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